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Abstract 


This  report  presents  the  computational  study  on  the 
Gaussian-beam  propagation  in  the  bio-optical  material 
bacteriorhodopsin  with  considerations  of  the  material’s 
intensity-dependent  absorption  and  refractive  index 
modulation.  The  beam  focusing  size,  focusing  position, 
intensity  change  in  the  material  and  their  dependence  on  the 
incident  beam  parameters  are  simulated. 


Iv 


I. 


Introduction 


Three-dimensional  (3-D)  optical  memories  [1-3]  show  great  promise  for  data  storage  and 
retrieval  because  of  its  high  capacity.  It  is  regarded  as  a  potential  candidate  for  future  high- 
performance  digital  applications,  such  as  multimedia  and  video  on  demand,  and  military 
and  industrial  object  recognition  applications.  Optical  memory  technology  poses  challenges 
for  scientists  to  search  for  an  ideal  photoactive  material.  The  photochromic  protein 
bacteriorhodopsin  (BR)  [4]  is  a  material  that  demonstrates  real-time  information  recording 
capabilities  and  has  a  promising  information-storage  potential  [5].  BR  undergoes  a 
complex  photocycle  on  the  absorption  of  light  energy.  The  two  relevant  photochemical 
states  are  called  bR  and  M  states  with  bR  being  the  ground  and  M  being  the  thermal 
intermediate  states,  respectively.  The  bR  state  has  a  strong  absorption  in  the  green-yellow 
region  of  the  visible  spectrum  (Xmax  =  568  nm).  The  M  state  exhibits  a  strong  blue-shifted 
absorption  spectrum  (Xmax  =  412  nm).  When  a  beam  is  traveling  in  this  material,  the 
absorption  and  phase  modulation  are  both  wavelength-  and  intensity-dependent.  Further 
propagation  of  the  beam  is  affected  by  these  light-induced  processes. 

In  both  two-photon  absorption  [7,8]  or  Fourier  transform  hologram  types  of  3-D 
optical  memories,  the  beam  focusing  property  plays  an  essential  role  in  determining  the 
information  storage  density  and  addressing  accuracy.  To  be  specific,  the  storage  density  in 
BR  or  other  storage  materials  depends  on  how  small  a  light  spot  can  be  formed  in  the 
material.  The  smaller  the  spot,  the  higher  the  possible  storage  density.  In  addition,  the 
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location  of  the  focused  light  spot  needs  to  be  accurately  determined  for  optimum 
distribution  and  precise  addressing  of  stored  data.  Because  of  the  nonlinear  wavefront 
modulation  of  BR,  the  above  issues  become  quite  involved  and  depend  on  the  wavelength, 
beam  power  and  divergent  angle,  as  well  as  the  material  properties. 

In  this  report,  we  present  the  computational  study  on  the  steady  state  behavior  of  a 
Gaussian  beam  as  it  passes  through  a  BR  film.  We  find  that  the  beam  waist  and  its  location 
in  BR  are  all  greatly  changed  due  to  the  nonlinearity  of  BR.  In  addition,  the  maximum 
intensity  no  longer  takes  place  at  the  beam  waist.  All  these  properties  may  be  important 
when  the  material  is  used  for  volume  Fourier  transform  holographic  memoiy  and/or  two- 
photon  absorption  optical  memory.  The  study  of  Gaussian  beam  propagation  is  useful  for 
future  studies  on  BR  as  an  optical  memory  medium.  This  study  can  also  be  used  as  a 
reference  to  other  absorption  and  index  modulation  recording  materials,  where  the  recording 
density  is  also  ruled  by  the  resolution  the  imaging  system  can  achieve.  In  previous  reports, 
[6,9]  we  have  presented  the  experimental  determination  of  the  effective  nonlinearity  of  BR 
film.  The  parameters  for  the  absorption  of  the  BR  film  were  also  reported.  [10]  Our 
simulations  are  based  on  these  parameters. 


II.  Optical  properties  of  bacteriorhodopsin  material 

Bacteriorhodopsin  is  a  light  transducing  protein  found  in  the  purple  membrane  of  the 
organism  Halobacterium  salinarium.  This  bacterium  grows  in  salt  marshes,  where  the 
concentration  of  salt  is  six  times  that  of  seawater.  A  BR  molecule  pumps  a  proton  across  the 
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bacterium’s  membrane  when  it  absorbs  light,  thus  generates  a  chemical  and  osmotic 
potential  that  serves  as  an  alternative  source  of  energy. 

The  number  of  times  that  BR  can  be  cycled  photochemically  -  switched  optically 
from  one  state  to  another  -  is  about  ten  million  times,  a  value  considerably  higher  than  that 
of  most  synthetic  photochromic  materials,  which  can  cycle  form  100  to  10,000.  Therefore, 
the  common  perception  that  biological  materials  are  too  fragile  for  use  in  optical  devices 
does  not  apply  to  BR. 

The  absorption  of  light  energy  by  theretinyl  chromophore  initiates  a  complex 
photocycle  characterized  by  several  spectroscopically  distinct  thermal  intermediates  and  a 
cycle  time  of  10  ms.  The  large  shift  in  absorption  maxima  accompanying  the  BR  to  M 
phototransformation  generates  a  corresponding  linear  change  in  the  refractive  index  which 
is  key  to  many  optical  processing  applications. 

BR  has  been  reported  to  be  suitable  for  three-dimensional  (3-D)  optical  memories 
[8].  Two-photon  absorption  is  regarded  the  one  of  the  prime  factors  for  the  information 
recording  in  3-D  optical  memories.  The  forward  reaction  from  bR  state  to  M  state  is  driven 
by  two  1 , 140-nm  photons,  and  the  reverse  reaction  from  M  state  to  bR  state  is  driven  by  two 
820-nm  photons. 
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Concept  of  a  3-D  optical  memory 


Figure  1.  A  3-D  optical  memory  information  storage  and 
retrieval. 


Figure  1  shows  a  typical  system  for  two-photon  absorption  3-D  optical  memory.  The 
object  beam  is  used  as  the  information  carrier  when  it  transmits  through  the  transmission 
input  spatial  light  modulator  (SLM).  The  to-be-recorded  information  is  electrically 
written  on  the  SLM.  The  object  is  imaged  onto  the  recording  medium  by  a  combination 
of  lenses  including  a  variable  focus-length  lens  on  the  Fourier  plane.  The  reference  beam 
is  the  other  arm  of  the  system  for  the  two-photon  absorption  photochromic  recording  of 
information.  Reference  beam  is  also  controllable  by  a  combined  lens  consisting  of 
conventional  lenses  and  an  electro-optic  lens.  The  reference  beam  may  produce  a  sheet  of 
light  to  cover  the  image  plane  of  the  SLM.  Two-photon  absorption  occurs  at  the 
interaction  section  of  the  two  beams.  Thus  the  two-dimensional  information  is  recorded  at 
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once.  When  the  variable  focus-length  lens  adjusts  its  focal  length,  the  image  plane  of 
SLM  shifts  along  the  optical  axis  in  the  recording  medium.  The  light  sheet  traces  the 
image  plane  through  a  scanning  device.  Thus  information  can  be  recorded  in  the  medium 
page  by  page.  The  page  information  storage  capacity  in  the  media  is  limited  by  the 
thickness  of  the  sheet-of-light. 

Since  there  is  a  physical  limitation,  although  the  process  is  massive  parallel,  the 
sheet  of  light  concept  does  permit  a  high  recording  density.  It  is  impossible  to  produce  a 
sheet  of  light  with  thin  thickness.  But  it  is  easy  to  achieve  a  fine  line  of  light.  Thus  as  an 
alternative,  we  proposed  the  concept  of  scanning  line-of-light.  The  information  recording 
is  performed  line  by  line  instead  of  page-by-page.  Although  parallelism  is  sacrificed 
compared  with  the  later,  the  line-by-line  recording  still  maintains  a  very  high  parallelism 
as  compared  with  the  current  serial  method  by  today’s  computer  industry.  And  because 
the  line-of-light  can  be  very  thin  (only  limited  by  the  real  system  design),  high  recording 
density  is  surely  promised. 

Three-dimensional  memory  can  also  be  implemented  by  using  focusing  point  of 
light  in  both  of  reference  and  object  beams.  This  will  give  a  memory  using  serial 
recording  and  information  retrieval.  We  can  achieve  the  highest  recording  density  but 
sacrifice  operation  time. 

Three-dimensional  memories  are  also  reported  with  other  techniques,  such  as 
wavelength  and  angular  multiplexing.  The  storage  space  for  wavelength  and  angular 
multiplexing  are  overlapped  in  the  recording  medium,  thus  suffering  low  recording  and 
diffraction  efficiencies.  Whichever  scheme  is  chosen  for  the  3-D  optical  memory 
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implementation,  the  laser  beam’s  propagation  is  the  essential.  Thus  a  study  of  the 
properties  of  Gaussian  beam  propagation  in  the  concerned  recording  medium  is  of  great 
interest. 


IV.  Modeling  and  Analysis 

A  Gaussian  beam  in  a  uniform  linear  medium  is  characterized  by  a  Gaussian  amplitude 
profile  and  a  quadratic  phase  distribution,  which  is  mathematically  expressed  as  [1 1] 


E{x,y,z)  =  E,-:^  exp  - 


W^iz) 


exp  ikz  +  ik 


2R(z) 


-ifpiz) 


where  k  is  the  wave  number  of  the  Gaussian  beam  in  the  medium.  Wo  is  the  waist  size,  and 
Eo  is  the  field  amplitude  at  the  waist  position. 


W\z)  =  W,^ 


(2) 


_  _KW,\  kW,\ 

X  r~ 


(3) 


R{z)=  z 


1  + 


(4) 


^(z)  =  tan‘'  —  ,  (5) 

no  is  the  intrinsic  index  of  refraction  of  the  medium.  In  a  nonlinear  medium,  however,  the 
amplitude,  phase,  and  divergence  angle  will  be  modulated  by  the  light  field  itself. 

BR  is  a  material  that  demonstiates  saturable  absorption  and  phase  modulation 
properties.  The  intensity-dependent  exponential  absorption  constant  at  the  steady  state  can 
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be  expressed  by  [10] 


a(/)=  A^a, 


l  +  2a^Tllhv 
l  +  (cT,  +a^)Tl/hv 


=  oc^- 


gl 


°  1  +  //// 


(6) 


where  =  Na^ ,  g  =  /Va,((T, -a,)vlhv  ,  =  hv/{o^  is  the  saturation  intensity, 

N  is  the  active  density  of  BR  molecules,  (7/  and  02  are  the  absorption  cross  sections  for  the 
nonradiative  transitions  bR=^  M  and  M  =>  fc/? ,  T  is  the  relaxation  time  for  the  transition 


M  =^bR,  vh  the  frequency.  The  saturation  intensity,  Is,  depends  on  the  lifetime  of  the  M 
state  and  the  absorption  cross-sections  of  the  BR  material.  For  the  laser  wavelength  of  632.8 
nm,  the  typical  parameters  for  BR-doped  polymer  films  in  the  above  equation  are  [9,10] 

Gq  =  8.2  cm”' 

I,  =  4.1  mW/cm'^  -  O) 

g  =  0.83  cmlmW 

For  the  laser  wavelength  of  5 14.5  nm,  the  typical  parameters  are  [9,10] 

=  17.2  cm”' 

I,  =9.5  mWicm-.  (8) 

g  =1.23  cmImW 


For  BR  prepared  with  different  options,  the  parameters  may  have  values  different  from  Eqs. 
(7)  and  (8).  The  saturation  intensity,  for  example,  is  a  function  of  the  lifetime  of  the  M  state, 
which  can  be  changed  by  chemical  enhancement  or  other  method  [6].  However,  the 
simulations  with  these  numbers  yield  typical  BR  behavior.  The  intensity-induced  phase 
change  under  low  intensities  can  be  expressed  by  [9,12] 


n  =  /Iq  +  An(/)  =  no 

where /io=  1.5,  =-1.25x10”’  cm^/mW  for  632.8  nm  wavelength  and  =+'^•25x10”’ 
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cm7mW  for  514.5  nm  wavelength,  n?  may  have  different  values  depending  on  the 
individual  preparation  process  of  the  BR  film  [9].  Eq.  (9)  may  be  held  correct  when  I  is 
small.  When  it  is  getting  larger,  the  index  change  will  reach  to  nonlinear  maximum  until  the 
thermal  effect  take  place.  [10]  At  that  point  we  may  assume  a  formula  for  the  relationship 
between  the  index  change  and  the  intensity  [9, 14] 

(10) 

1  +  — 

h 

The  maximum  of  An  is  assumed  [13,14]  to  be  -l.OxlO’’  for  red  light  and  l.OxlO"'’  for 
green  light,  respectively.  Both  and  are  functions  of  molar  refraction  coefficients  and 
concentrations  of  the  M  and  bR  states.  Therefore,  with  and  the  maximum  An  known, 
we  have  the  estimation  of 


=  8.0x10-'  mW/cm'. 


(11) 


The  propagation  of  a  coherent  paraxial  beam  through  a  distance  in  a  uniform  and 
linear  medium  can  be  obtained  by  using  the  spatial  transfer  function  for  propagation  of  d 
distance  in  that  medium.  We  modify  this  method  to  study  the  light  field  evolution  of  a 
Gaussian  beam  during  its  propagation  in  the  BR  film.  We  treat  the  propagation  in  two  steps: 
first,  the  nonlinear  amplitude  and  phase  modulation  due  to  the  BR  material  and  second,  the 
diffraction  effect  of  wave  propagation.  We  assume  that  the  intensity  distribution  of  the 
Gaussian  beam  in  the  input  plane  of  a  differential  layer  of  thickness  d  determines  the 
absorption  and  phase  modulation  properties  of  this  layer.  So  the  nonlinear  wave 
modification  by  this  layer  can  be  modeled  as  an  infinitely  thin  layer  of  amplitude 
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transmittance  t,  that  is  expressed  as 

/(jc,>';z)  =  exp|-^a(|£(x,y;z)|')/  +  /^oAnj£(x,3;;z]|^)/|,  (12) 

where  a  constant  phase  factor  has  been  ignored.  The  electrical  field  has  absorbed  a  consttint 
factor  so  that  we  denote  I  =\Ef .  Taking  Eq.  (1)  as  the  boundary  condition,  the  effective 
light  field  that  will  propagate  through  the  layer  of  thickness  d  is  given  by 

Eeff  (•*’  y'^^)=  y\z)E{x,y,z).  (13) 

This  modified  wavefront  then  propagates  through  the  distance  r/  as  if  it  was  a  “uniform 
linear  medium”  whose  refractive  index  is  constant.  In  other  words,  only  diffraction  effects 
are  considered  for  E^x,y\z). 

A  fast  Fourier  transform  is  used  to  convert  the  effective  spatial  light  distribution, 
Eej^x,  y;  z),  into  the  angular  plane-wave  spectrum,  %/p,  q\  z),  which  satisfies 


'E^{p.q\z)=  j  \E^ff{x,y\z)txp[-i{px-\-qy)]dxdy 

—00—00 
oo  oo 

E,ff{p,q\z)=  J  l'P,ff{p,q;z)Qxp[i{px  +  qy)]dpdq 


(14) 


The  spectrum  distribution  is  multiplied  by  the  spatial  transfer  function  of  propagation  of 
distance  d.  The  spatial  transfer  function  of  this  paraxial  propagation  for  this  distance  can  be 
expressed  as 

H{p,q,d)=&xp' 

=  exp- 

where  p  and  q  are  coordinates  for  spatial  angular  frequencies  and  Xq  is  the  fioe-space 


2k  J 

\iXod(p^+q^)] 

I  4;!n„ 


(15) 
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wavelength  of  the  incident  beam.  A  followed  inverse  Fourier  transform  gives  us  the  light 
distribution  after  the  ^/-distance  propagation.  Thus 

E{x,y,z+d)=F  ''[H{p,q-,d)F  y; z)]} 

=  F  ''{H{p,q-,d)F  [r(x,y;z)£(j:,y;z)]}, 

where  F  and  F  '*  denote  Fourier  and  inverse  Fourier  transforms,  respectively.  The  near- 
field  Fresnel  diffraction  is  obtained  by  repeating  the  above  process  layer  by  layer  to  find  out 
the  wave  evolution.  The  separated  treatments  have  little  effect  on  the  accuracy  of  the  field 
expression  when  d  is  small  enough.  This  assumption  is  tested  in  the  simulation  by  the  fact 
that  changing  d  in  a  certain  range  does  not  affect  the  final  field  distribution. 

The  focusing  property  of  the  beam  in  the  medium  is  an  important  parameter.  It 
determines  the  information  storage  density  in  that  medium.  We  started  with  a  waist  size  of  1 
pm,  a  practical  size  for  any  meaningful  optical  storage.  Because  of  the  large  divergent  angle 
for  this  waist,  we  simulated  a  one-dimensional  array  in  order  to  cope  with  the  huge 
computing  task  without  compromising  the  generality.  For  one-dimensional  case,  the 
complex  amplitude  of  the  electrical  field,  expressed  by  Eq.  (1)  can  be  replaced  by 


E{x\z)=E^ 


'w(z) 


exp<  ikz  -  X* 


W^{z)  2R{z) 


the  parameters  are  the  same  as  defined  in  Eqs.  (2)-(5).  The  propagation  factor  now  has  p  as 
its  sole  variable.  In  the  simulation,  we  use  laser  beam  power,  denoted  by  PWR,  as  a 


PWR  2 


parameter.  The  one-dimensional  Gaussian  beam’s  intensity  will  be/^  = -  —  while  it 

V;r 


2'x.PWR 

*•’’  ^0  ~ . -  for  a  two-dimensional  case.  A  one-dimensional  Gaussian  beam  is 
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generated  and  put  into  a  4096-point  one-dimensional  array.  The  BR  film  is  assumed  with  a 
thickness  of  600  p,m.  The  film  is  divided  into  300  thin  layers.  We  assume  that  the  waist  size 
is  focused  onto  a  position  360  [im  deep  into  the  BR  if  no  absorption  and  phase-change  are 
present. 


V.  Computer  simulations 

Simulations  were  performed  with  two  typical  laser  wavelengths.  One  is  He-Ne’s  632.8  nm, 
the  other  is  Ar  ion’s  514.5  nm.  They  present  different  refractive  index  change  properties. 
First,  we  consider  the  effect  of  the  nonlinearity  on  the  waist  size  of  a  Gaussian  beam  that 
propagates  through  a  BR  film.  Besides  determining  the  information  density  for  bitwise 
storage,  this  is  also  useful  for  image  hologram  where  a  high  density  recording  requires  a 
small  recording  area  onto  the  BR  film.  During  this  simulation,  we  choose  the  power  of  the 
laser  beam  to  be  1.0,  0.1,  0.01  and  0.001  mW  and  keep  the  saturation  intensity,  /„ 
unchanged.  The  incident  waist  position  was  focused  360  p.m  deep  into  the  BR  film  provided 
there  would  be  no  nonlinearity  presented.  Because  this  is  a  one-dimensional  simulation,  we, 
when  using  intensity  or  power  of  the  beam,  assume  that  the  size  of  the  beam  in  the  other 
dimension  is  1  mm. 

The  relative  change  of  beam  waist  size  in  BR  as  a  function  of  the  incident  beam 
waist  size  at  a  variety  of  laser  beam  powers  is  shown  in  Figure  2.  When  the  incident 
Gaussian  beam  is  with  a  low  power  (0.01  mW  and  0.001  mW),  the  actual  waist  size  in  BR 
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Figure  2.  Influence  of  nonlinearity  on  the  waist  size  of  a 
Gaussian  beam  for  (a)  red  iight  of  632.8  nm  and  (b)  green  light 
of  514.5  nm.  The  relative  change  is  defined  as  the  actual  waist 
size  minus  the  original  incident  waist  size. 
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is  almost  the  same  as  the  incident  beam  size.  The  actual  focused  size  is  always  smaller  than 
that  of  the  incident  when  the  incident  waist  is  large  (i.e.  incident  waist  size  is  about  10  |xm 
or  above).  This  is  because  the  phase  change  effect  is  smaller  than  that  of  the  amplitude 
absorption.  The  intensity-dependent  absorption  makes  the  waist  size  smaller  because  low 
intensity  suffers  a  large  absorption.  When  the  beam  power  is  high  enough  (0.1  mW  and  1.0 
mW),  the  actual  waist  differs  greatly  from  the  incident  beam  waist  size.  For  a  very  small 
waist  (1  or  2  |im)  the  intensity  of  the  beam  varies  dramatically  during  its  propagation 
because  of  its  large  divergent/convergent  angle.  Thus  the  major  nonlinear  effect  affecting 
the  waist  size  can  only  happen  in  a  very  limited  space  near  the  beam  waist.  Therefore  the 
effect  on  the  waist  size  change  is  small.  But  if  waist  size  is  large  enough,  the  diverging 
angle  is  not  large  and  thus  the  nonlinear  effect  can  happen  in  a  relatively  long  distance  of 
beam  propagation.  The  negative  index  change  for  He-Ne  laser  at  a  high  intensity  makes  the 
converging  beams  divergent  and  reaches  a  waist  phase  distribution  (i.e.  uniform  phase 
distribution)  earlier  than  expected.  Thus  a  larger-than-expected  wedst  size  yields.  Even  in 
the  low  incident  waist  regions  of  0.01  mW  or  0.001  mW  this  diverging  effect  is  noticeable. 
The  positive  index  change  for  514.5  nm  light  enables  the  green  light  to  be  focused  in  a  way 
opposite  to  that  of  the  red  light  and  forms  a  waist  size  much  smaller  than  that  of  the  incident 
wave.  The  most  interesting  results  are  with  the  power  of  1.0  mW.  As  shown  in  Figure  2(b), 
the  relative  change  of  the  waist  size  is  almost  proportional  to  the  incident  beam  waist  size. 
As  the  incident  waist  size  increases,  the  actual  waist  size  in  BR  does  not  increase  that  much 
and  is  kept  around  just  a  few  microns.  It  is  noted  that  when  the  incident  waist  size  is  larger 
than  16  pm,  the  actually-focused  beam  size  increases  as  the  incident  waist  size  increases. 
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Figure  3.  Focusing  position  as  a  function  of  the  incident 
beam's  waist  size  for  (a)  red  light  of  632.8  nm  and  (b)  green 
light  of  514.5  nm.  The  incident  light  is  supposed  to  be  focused 
at  360  (xm  if  no  nonlinear  absorption  and  phase  modulation  are 
present  in  BR. 


This  Phenomena  is  brought  out  because  we  only  consider  a  film  of  600  jxm.  The  actual 
focus  position  is  out  of  the  film  and  our  data  just  includes  the  minimum  value  inside  film. 
This  can  be  clearly  seen  with  the  help  of  Figure  3(b). 

The  relationship  between  the  focusing  position  and  the  incident  waist  size  was  also 
simulated,  as  shown  in  Figure  3.  As  the  plots  show,  the  actual  beam  waist  location  in  BR 
can  differ  greatly  from  that  of  the  original  incident  beam,  and  the  differences  depend  on  the 
incident  wavelength,  incident  waist  size,  and  power.  For  the  red  light  of  632.8-nm 
wavelength,  the  actual  focusing  position  can  be  either  ahead  of  or  beyond  the  incident  waist 
position.  For  a  low  beam  power  (0.01  mW  and  0.001  mW),  the  beam  intensity  is  low  and 
the  intensity-induced  phase  modulation  slightly  reduces  the  converging  beam’s  convergence 
but  not  strongly  enough  to  make  the  beam  divergent  quickly.  Thus  actual  focusing  position 
is  moved  a  little  bit  towards  the  front  surface  of  the  BR  film.  Meanwhile  the  absorption 
modulation  dominantly  affects  the  focusing  position.  The  weaker  peripheral  part  of  the 
Gaussian  beam  is  strongly  absorbed  as  the  beam  propagates  and  this  effect  always  moves 
the  position  away  behind  the  incident  waist  position.  When  the  power  is  high,  the  diverging 
nature  of  strong  index  modulation  makes  the  converging  laser  beam  divergent  before  it  can 
reach  the  incident  waist  position  and  yields  a  large  actual  focusing  beam  size.  Exceptions 
are  with  the  1  to  3  pm  region  of  the  curve  for  1.0  mW  where  the  beam  is  with  a  large 
converging/diverging  angle.  Thus  at  the  front  surface  of  the  BR  film  the  beam  size  is  large 
and  the  intensity  is  relatively  low.  The  intensity  distribution  from  the  front  surface  to  the 
focus  point  changes  dramatically.  First  the  intensity  is  low;  thus  the  absorption  takes  a  main 
role.  Then  the  intensity  becomes  high,  and  the  index  change  dominates  the  effect.  The 
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above-mentioned  exceptions  are  the  combination  with  the  absorption  taking  a  main  part. 
For  the  green  light  of  514.5  nm  as  shown  in  Figure  3(b),  most  of  the  focusing  position  for  a 
large  waist  stops  at  600  p,m  as  the  incident  waist  size  increases.  This  is  because  we  only 
consider  a  BR  film  that  is  600-|xm  thick.  The  focused  point  have  reached  out  the  back 
surface  of  the  BR  film,  hence  the  minimum  beam  size  within  the  BR  film  is  located  at  the 
back  surface  of  the  BR  film.  Where  the  curve  reaches  the  600-p.m  boundary  corresponds  to 
the  turning  point  in  Figure  2(b)  of  the  corresponding  curves.  When  the  power  is  low  (0.01 
mW  and  0.001  mW),  absorption  takes  the  main  role  influencing  the  focusing  position.  The 
curves  are  similar  to  what  Figure  3(a)  shows  for  the  He-Ne  low  beam  power.  When  the 
power  increases  (0.1  mW  to  1.0  mW),  the  converging  nature  of  the  modulation  makes  the 
Gaussian  beam  a  little  more  convergent  and  reaches  focused  position  earlier  than  the 
incident  waist  position.  The  increase  of  incident  beam  waist  size  reduces  the  index  phase 
modulation  and  strengthens  the  absorption  modulation  effect.  Therefore  the  focusing 
position  moves  back  to  the  back  surface  of  the  BR  film.  The  results  of  Figures  2  and  3 
indicate  that  the  beam  waist  size  and  position  will  change  in  BR  in  a  nonlinear  manner  and 
some  of  these  changes  are  so  big  that  they  should  be  considered  when  memory  position 
allocation  is  considered,  hi  other  words,  pre-recording  calibration  is  required  in  order  to 
record  reliable  data  in  the  BR  film. 

When  we  look  at  the  incident  waist  position,  that  is,  360  pm,  the  actual  beam  sizes 
are  different  from  that  of  the  incident  beam,  as  shown  in  Figure  4.  The  vertical  axis  denotes 
the  actual  beam  width.  For  the  He-Ne  wavelength  beam  width  is  always  larger  than  the 
incident  waist  size.  The  largest  error  occurs  when  the  waist  size  is  around  5  pm.  The  actual 
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Rgure4.  Actual  beam  width  at  the  original  incident  beam 
waist  position  for  (a)  wavelength  of  632.8  nm  and  (b) 
wavelength  of  514.5  nm. 
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beam  width  almost  double  that  of  the  incident  beam.  But  the  case  for  an  Ar  green  line  is 
more  complicated.  In  most  of  the  case,  the  beam  width  at  the  expected  focusing  position  has 
a  reduced  waist  size.  A  larger-than-expected  beam  size  occurs  only  when  the  expected  waist 
size  is  small  and  beam  power  is  high. 

Because  of  the  presence  of  both  absorption  and  phase  modulation  in  BR,  the  highest 
intensity  does  not  occur  at  the  position  of  the  beam  waist.  The  difference  between  these  two 
positions  as  a  function  of  incident  waist  size  is  depicted  in  Figure  5.  Because  this  difference 
is  caused  by  the  absorption,  the  highest  intensity  generally  occurs  before  the  waist  position. 
There  is  only  one  exception  in  the  graph  of  the  He-Ne  light  with  the  Gaussian  beam  power 
to  be  1 .0  mW  and  the  target  waist  is  4  [xm.  The  beam  intensity  profiles  for  this  exception  are 
drawn  in  Figure  6.  Attention  is  given  to  three  situations,  the  intensity  distribution  at  the 
front  surface  of  the  BR,  that  at  the  focusing  position  and  that  at  the  maximum  central 
intensity  location.  The  first  case  is  indicated  by  a  thick  line  in  the  figure  while  the  other  two 
are  represented  by  a  medium  line  and  a  thin  line,  respectively.  The  first  case  is  an  ideal 
Gaussian  distribution.  But  due  to  the  nonlinear  modulation,  the  beam  intensity  profiles  for 
the  second  and  third  cases  do  not  maintain  a  Gaussian  shape,  though  very  alike.  When 
judged  by  the  Gaussian  waist  size  evaluation  criteria,  i.e.  the  distance  between  the  locations 
where  the  intensity  drops  to  1/e^  of  the  central  intensity,  the  second  case  is  the  best.  There 
are  also  two  exceptions  for  the  green  light  at  10  (im  and  1 1  |xm.  These  are  merely  caused  by 
simulation  accuracy.  When  the  beam  power  is  high,  as  shown  in  Figure  5,  the  absorption 
effect  does  not  take  main  part  and  the  difference  is  not  large.  When  the  beam  power 
becomes  higher  and  higher,  the  difference  gets  larger  and  larger.  There  are  turning  points  on 
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Figures.  The  difference  between  the  beam  waist  position 
and  the  highest  intensity  position  for  (a)  632.8  nm  wavelength 
and  (b)  514.5  nm  wavelength. 
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Figure  6.  Profiles  of  the  beam's  intensity  distribution  for  the 
red  light  at  front  surface  of  the  BR  film  and  positions  inside  the 
BR  film  where  we  found  the  actual  focus  and  maximum  central 
intensity.  The  incident  beam  power  is  1.0  mW  and  the  waist 
size  is  4  urn. 

the  low  beam  power  curves.  For  the  He-Ne  light  beam,  turning  points  occur  at  a  10-|jjn 
incident  waist  size  with  beam  powers  of  O.l  mW  and  0.01  mW.  That  is  when  the 
maximum-central-intensity  point  reaches  the  front  surface.  This  position  can  not  be  moved 
any  further.  Therefore  further  increase  of  the  incident  beam  waist  size  see  a  slow  increase  of 
the  difference  between  the  focusing  position  and  the  maximum  central  intensity  position. 
The  situation  is  more  complicated  for  green  light.  For  the  high  beam  powers  (1.0  mW,  0.5 
mW  and  0.3  mW),  there  is  little  difference  between  the  locations  of  focusing  position  and 
maximum  central  intensity.  As  the  incident  beam  waist  increases,  both  these  positions  move 
towards  the  back  surface  of  the  BR  film.  For  the  0. 1  mW  beam  power  curve,  the  increase  of 
the  incident  waist  size  moves  the  maximum  central  intensity  position  towards  the  back 
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surface  of  the  BR.  But  when  incident  waist  is  large  than  1 1  |J.m,  the  absorption  effect  is  so 
strong  that  the  central  intensity  at  the  front  surface  of  the  BR  film  is  larger  than  any  central 
intensity  in  the  film,  and  thus  a  huge  jump  is  presented  on  this  curve.  To  present  this 
situation  clearly.  Figure  7  shows  the  beam  central  intensity  as  a  function  of  the  BR  film’s 
depth  for  a  number  of  incident  waist  sizes. 


Figure  7.  Central  intensity  via  BR  film  depth  relationship  for 
a  number  of  incident  beam  waist  sizes  at  a  beam  power  of  0.1 
mW  of  green  light.  Curves  from  top  to  bottom  are 
corresponding  to  an  incident  waist  size  of  8,  9,  10,  11  and  12 
pm. 

The  curves  from  top  to  bottom  are  corresponding  to  an  incident  waist  size  of  8,  9, 
10,  1 1  and  12  pm.  When  incident  waist  size  varies  from  8  pm  to  10  pm,  the  maximum 
central  intensity  (peak  value)  moves  towards  the  back  surface  of  the  BR  film.  When  the 
incident  waist  size  becomes  11  pm,  although  the  peak  intensity  moves  in  the  same 
direction,  the  maximum  central  intensity  occurs  at  BR’s  front  surface.  The  fact  that  the 
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maximum  intensity  no  longer  occurs  at  the  beam  waist  is  significant  in  3-D  optical  memory 
applications.  In  two-photon  optical  memory,  the  maximum  intensity  location  is  where  the 
two-photon  absorption  happens.  In  Fourier  volume  holographic  memory,  the  location  of  the 
maximum  intensity  should  be  considered  for  optimum  interference  modulation. 

For  the  632.8-nm  wavelength  the  focused  beam  size  is  generally  larger  than  the 
incident  beam  waist  while  for  the  5l4.5-nm  wavelength  the  focused  beam  size  can  be 
smaller  than  the  waist  size  of  the  incident  beam.  This  property  for  the  green  light  may  be 
used  to  achieve  a  beam  size  smaller  than  the  diffraction  limit.  We  simulated  the  focused 
beam  size  as  a  function  of  the  illuminating  beam  power.  The  simulation  results  for  an 
incident  waist  position  of  100  |xm  and  360  pm  are  shown  in  Figures  8  and  9,  respectively. 
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Figure  8.  Waist  size  via  the  laser  beam  power  for  632.8  nm 
and  5 14.5  nm  for  incident  waist  sizes  of  3,  5,  7  and  9  pm.  The 
incident  waist  position  is  360  p,m  into  the  BR  film. 
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Figure9.  Waist  size  via  the  laser  beam  power  for  632.8  nm 
and  514.5  nm  for  incident  waist  sizes  of  3,  5,  7  and  9  |i.m.  The 
incident  waist  position  is  100  )xm  into  the  BR  film. 


Both  red  and  green  wavelengths  are  considered  with  incident  waist  sizes  of  3,  5  and 
7  pm.  The  three  top  lines  are  for  red  light  and  the  three  bottom  lines  are  for  green  light.  For 
both  the  green  and  red  light,  there  exists  a  beam  power  at  which  an  extremum  focus  size  is 
formed.  The  positive  index  change  of  the  green  light  provides  a  converging  lensing  effect 
and  increases  the  converging  angle  of  the  Gaussian  beam,  hence  a  smaller  beam  size  is 
formed.  The  increase  of  the  beam  power  intensifies  the  lensing  effect  and  reduces  the 
focused  beam  size.  For  the  red  light  there  is  a  maximum  focused  size.  As  the  beam  power 
increases  from  null,  the  increased  index  change  makes  the  red  to  have  larger  and  larger 
focused  size  and  the  green  light  smaller  and  smaller  one,  as  those  curves  shown  in  Figure  2. 
But  when  the  power  is  even  higher,  the  saturation  in  absorption  and  index  change  expressed 
by  Eqs.  (6)  and  (10)  reduces  the  modulation  on  the  incident  beam  and  causes  the  red  light’s 


23 


focused  size  to  decrease  and  the  green  light's  focus  size  to  increase.  The  beam  power  to 
reach  an  extremum  waist  size  depends  on  the  incident  waist  size  and  the  incident  focusing 
position.  A  small  incident  waist  size  requires  a  smaller  beam  power  to  result  an  extremum 
of  focused  beam  size  than  that  for  a  large  incident  waist  size;  and  a  short  incident  waist 
position,  by  comparing  Figure  9  with  Figure  8,  can  also  provide  the  similar  result. 


VI.  Conclusions 

We  have,  with  the  632.8-nm  He-Ne  laser  wavelength  and  5l4.5-nm  Ar  wavelength, 
numerically  investigated  the  Gaussian  beam  propagation  property  through  a  thick  BR 
film.  It  is  shown  that  the  focusing  position  and  the  waist  size  of  the  Gaussian  beam  are 
both  affected  by  the  nonlinearity  of  the  BR  material.  Thus  when  using  BR  as  the  optical 
memory  media  these  parameters  should  be  taken  into  account  in  order  to  get  an  optimum 
and  reliable  recording.  The  method  can  also  be  applied  to  analyze  other  intensity- 
dependent  absorption  and  phase  modulation  materials.  From  this  simulation,  it  is  also 
concluded  that  for  a  very  large  converging  angle,  i.e.  an  incident  waist  size  under  1  |im, 
the  focused  beam  size  and  position  are  not  greatly  influenced  by  the  nonlinearity  of  the 
material. 
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